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We analyzed the prokaryotic commurnity structure of a saltern pond with 21 % total salts 
located in Isia Cristina, Huelva, Southwest Spain, close to the Atlantic ocean coast. For this 
purpose, we constructed a metagenome (designated as IC21 ) obtained by pyrosequencing 
consisting of 486 Mb with an average read length of 397 bp and compared it with other 
metagenomic datasets obtained from ponds with 19, 33, and 37% total salts acquired 
from Santa Pola marine saltern, located in Alicante, East Spain, on the Mediterranean 
coast. Although the salinity in IC21 is closer to the pond with 19% total salts from Santa 
Pola saltern (designated as SSI 9), IC21 is more similar at higher taxonomic levels to the 
pond with 33% total salts from Santa Pola saltern (designated as SS33), since both are 
predominated by the phylum Euryarchaeota. However, there are significant differences at 
lower taxonomic levels where most sequences were related to the genus Halorubrum in 
IC21 and to Haloquadratum in SS33. Within the Bacteroidetes, the genus Psychroflexus 
is the most abundant in IC21 while Salinibacter dominates in SS33. Sequences related to 
bacteriorhodopsins and halorhodopsins correlate with the abundance of Haloquadratum 
in Santa Pola SSI 9 to SS33 and of Halorubrum in Isia Cristina IC21 dataset, respectively. 
Differences in composition might be attributed to local ecological conditions since IC21 
showed a decrease in the number of sequences related to the synthesis of compatible 
solutes and in the utilization of phosphonate. 
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INTRODUCTION 

Hypersaline habitats are characterized by high sah concentra- 
tions, in addition to other features, such as high or low tempera- 
tures, high pH, and/or low oxygen concentrations (lavor, 1989; 
Rodriguez- Valera, 1993). Hypersaline environments are often 
aquatic systems (thalassohaline, of marine origin, or athalassoha- 
line, formed by dissolution of mineral salt deposits of continental 
origin) or saline soils (Walsh et al., 2005; Ventosa, 2006; Ventosa 
et al., 2008), but they are also represented by salt deposits, some 
desert plants, oilfield brines and a variety of salted foods, from 
seasoned fish or meat to fermented foods as well as animal hides 
(Grant et al, 1998; Ventosa, 2006). The best studied hypersaline 
habitats are aquatic hypersaline systems, such as salt lakes and 
salterns. 

Salterns are excellent models for studying the ecology and 
diversity of microorganisms, given that they are composed by a 
series of ponds with widely different salinities that concentrate 
salt from seawater to the point of saturation and precipitation. 
Most saltern studies have been performed on the saturated brine 
crystallizer ponds (Anton et al, 1999, 2000; Benlloch et al., 2001; 
Pasic et al, 2005; Pasic et al, 2007; Legauh et al., 2006; Oh et al, 
2010). Some comparative reports performed on crystaUizers 
ponds, as Pasic et al. (2007), compared haloarchaeal communities 



from two Adriatic solar salterns, showing differences in the 
microbiota and that climate could play a role in the microbial 
community structure. Oh et al. (2010) examinated the diversity 
of Haloquadratum and other haloarchaea in three coastal, but 
geographically distant saltern crystallizer ponds in Australia. 
The great majority of the 16S rRNA gene sequences recovered 
from these crystaUizers were related to H. walsbyi and diverged 
by less than 2% from each other, and from the type strain 
of this genus (strain C23). However, our knowledge about 
intermediate salinity ponds is limited. Benlloch et al. (2002) 
analyzed 16S rRNA sequences by DGGE of three salt ponds (8, 
22, and 32% total salts, respectively) from Santa Pola saltern 
in Eastern Spain. Most bacterial sequences in the 8% salt 
pond were related to organisms of marine origin belonging to 
representatives of the classes Alpha-, Beta-, Gamma-, and Epsilon- 
Proteobacteria, and the phyla Bacteroidetes, Actinobacteria, and 
Cyanobacteria. In the 22% salt pond were found Alpha- and 
Gamma-Proteobacteria, Cyanobacteria, and Bacteroidetes, and 
most of them were related to specialized halophiles. From 
the 32% salt pond, the only Bacteria found were sequences 
that clustered with Salinibacter ruber, an extremely halophilic 
Bacteroidetes. And in those three different salinity ponds, most 
of the clones were related to cultured strains of the archaeal 
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class Halobacteria. A metagenomic study of an intermediate 
salinity pond from Santa Pola saltern (SS19), revealed a low 
presence of bacterial species of Halomonas, Chromohalobacter, 
or Salinivibrio (Ghai et al., 2011), which are commonly isolated 
from those habitats (Arahal et al, 2001; Arenas et al, 2009) 
but abundant metagenomic reads affiliated to Haloquadratum 
walsbyi and Salinibacter ruber were also found. A novel repre- 
sentative of nanohaloarchaeota {"Candidatus Haloredivivus") 
which was also found to be abundant, and its genome was 
partially assembled. The most abundant bacterium found in this 
19% salt pond appeared to be a gammaproteobacterium closely 
related to Alkalilimnicola and Nitrococcus. This microbe has been 
recently cultured and its genome sequenced (Leon et al., 2013; 
Lopez-Perez et al., 2013). Besides, Ghai et al. (2011) found a 
large number of sequences related to presumably non-halophilic 
bacterial genera and, a group of low G-l-C Actinobacteria typical 
to freshwater habitats. A recent metagenomic study on a pond 
with 13% salts of Santa Pola saltern showed a large microbial 
diversity representing seven major higher taxa: Euryarchaeota, 
Gammaproteobacteria, Alphaproteobacteria, Actinobacteria, 
Bacteroidetes, Verrucomicrobia, and Betaproteobacteria 
(Fernandez et al., 2014a). Community analysis of an inter- 
mediate salinity pond (18% total salts) from a saltern in Guerrero 
Negro, Mexico, demonstrated that the archaeal community was 
dominated by a single uncultured 16S rRNA phylotype with 
99% similarity to sequences recovered from a Tunisian saltern 
and the most abundant bacterial sequences were 99% similar 
to an uncultured gammaproteobacterial clone from the Salton 
Sea (Dillon et al, 2013). Boujelben et al. (2012) explored the 
prokaryotic community in several ponds from Sfax saltern in 
Tunisia. They showed that some phylotypes, such as those related 
to Haloquadratum or representatives of Bacteroidetes, displayed 
a strong dependence of salinity and/or magnesium concentra- 
tions and that temperature was a strong factor structuring the 
prokaryotic community in the pond with 20% salinity, but not 
in the crystallizer pond, due to the seasonal changes. However, 
a survey about six salt lakes in Inner Mongolia, China, and a 
salt lake in Argentina showed that archaeal biogeography was 
influenced by Na+, COj^, HCO^^, pH and temperature, and 
bacterial biogeography was influenced by Na+, Mg^+, HCO^", 
and pH as well as geographic distance (Pagaling et al, 2009; 
Grant et al, 2011). 



Therefore, in order to learn more about the community struc- 
ture of intermediate salinity ponds, we explored the phyloge- 
netic and taxonomic differences as well as metabolic profiles of 
two geographically distant habitats from Spain: the Isla Cristina 
saltern that gets its water from the Atlantic Ocean and the Santa 
Pola saltern on the Mediterranean Sea. 

MATERIALS AND METHODS 

SAMPLE COLLECTION, DNA EXTRACTION, AND SEQUENCING 

The metagenomic datasets analyzed in this study are derived 
from different saline systems: one dataset was obtained from a 
pond at the Isla Cristina saltern located on the Atlantic Ocean, 
in southwestern Spain with a salinity of 21% (Fernandez et al., 
2014b); four datasets were from the Santa Pola saltern located 
in eastern Spain on the Mediterranean Sea, of which three were 
from concentrator ponds (SS13, SS19, and SS33) and one from 
a crystallizer pond (SS37) with salt concentrations of 13, 19, 33, 
and 37%, respectively (Ghai et al, 2011; Fernandez et al., 2013, 
2014a); and two marine datasets, deep chlorophyll maximum 
from Mediterranean Sea (DCM3) and Mar Menor coastal lagoon 
(MM5) with salinities of 3.8 and 5%, respectively (Ghai et al., 
2010, 2012). AH databases were obtained using the same DNA 
extraction method and the samples were sequenced by pyrose- 
quencing 454 (Martin-Cuadrado et al, 2007; Ghai et al, 2011). 
The accession numbers for the deposited databases are shown in 
Table 1. 

COMPARATIVE ANALYSIS OF METAGENOMIC READS 

To estimate cumulative nucleotide differences between metage- 
nomic datasets, we carried out BLASTN searches of the complete 
set of sequences from every dataset vs. all the others. Bitscores of 
the top high-scoring segment pairs (HSPs) from every sequence 
from one set vs. another were summed to yield a cumulative 
pairwise bitscore value (CPBV) that was normalized and used to 
construct a distance matrix. CPBVs were normalized by dividing 
each one by the cumulative bistscore value derived from the 
BLASTN of one dataset vs. itself. The distance matrix was ana- 
lyzed using the Phylip package (Felsenstein, 1989) to obtain a 
neighbor-joining tree. 

G-l-C contents were computed using the program geecee in the 
EMBOSS package (Rice et al., 2000) and the amino acid frequency 
was calculated from a perl script. The metagenomic reads were 



Table 1 | Features of the different datasets from saline habitats used in this study. 



Datasets 


Salinity 

(%) 


Temperature 
(X) 


pH 


Number of 
reads 


Dataset 
size (Mb) 


Average read 
length (bp) 


Accession 
number/source 


References 




Deep Chlorophyll 


3.8 


15.9 


8.1 


1,204,321 


312 


259 


SRP002017 


Ghai et al., 2010 




Maximum (DCM3) 




















Mar Menor Coastal 


5 


19 


8.4 


730,997 


243 


335 


http://camera.calit2.net 


Ghai et al., 2012 




Lagoon (MM5) 




















Santa Pola Saltern (SS13) 


13 


29 


8.0 


1,443,593 


441 


305 


SRP028290 


Fernandez et al.. 


2013, 2014a 


Santa Pola Saltern (SS19) 


19 


30 


8.0 


1,315,302 


475 


361 


SRP007685 


Ghai et al., 2011 




Santa Pola Saltern (SS33) 


33 


30 


70 


842,872 


309 


367 


SRP028290 


Fernandez et al.. 


2013 


Santa Pola Saltern (SS37) 


37 


41 


71 


760,740 


309 


417 


SRP007685 


Ghai et al., 2011 




Isla Cristina Saltern (IC21) 


21 


25 


75 


1,223,923 


486 


397 


SRP029970 


Fernandez et al.. 


2014b 
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annotated using UniProtKB database released in December 2013 
(UniProt Consortium, 2014) through BLASTX search with a 
cutoff e-value of le-5. 

16S ribosomal RNA genes were identified by comparing the 
datasets against the RDP database version 11.1 (Cole etal.,2014). 
All reads that matched a 16S rRNA sequence with an alignment 
length of more than 100 bp and an e-value lower than le-5 against 
the database were extracted. The best hit that was not described 
as unknown or unidentified was considered a reasonable closest 
attempt for classifying the 16S rRNA sequences. Sequences were 
assigned to a specific genus if they shared >95% 16S rRNA 
sequence identity with a known representative. 

METAGENOMIC READS ASSEMBLY 

Assembly of the metagenomic reads with greater than 100 bp 
was performed using stringent criteria requiring an overlap of at 
least 80 bp, 99% identity and at most a single gap in the align- 
ment (using Geneious Pro 5.4). Next, assembled contigs that 
were less than 3 kb in length, and those with less than three 
predicted genes were discarded for further analysis. We retained 
only those contigs that provided consistent query hits to only 
single high level taxa (e.g., Alphaproteobacteria, Euryarchaeota, 
Bacteroidetes, Actinobacteria). To test if the assembly strategy pro- 
duced authentic contigs from known organisms, we manually 
identified all contigs that belonged to H. walsbyi, one of the abun- 
dant organism in the datasets. The criterion was that all genes 
from a putative H. walsbyi contig must return best hits from that 
genome. 

Tetranucleotide frequencies of the assembled contigs were 
computed using the wordfreq program in the EMBOSS package 
(Rice et al., 2000), and principal component analysis (PCA) was 
performed using the R package FactoMineR (Le et al, 2008). 

CONSTRUCTION OF PHYLOGENETIC TREES 

Maximum likelihood reference trees were constructed using 
RaxML as implemented in ARB software package (Ludwig et al., 
2004) using reference 16S rRNA gene sequences with near 
full length (>1300 nt) from cultured isolates. Later, partial 
16S rRNA gene assembled metagenomic sequences and closely 
related environmental uncultured 16S rRNA gene sequences were 
inserted into reference trees without altering tree topology using 
maximum parsimony criterion and a 50% base frequency fil- 
ter. Bootstrap values greater than 50% are indicated above nodes 
and the scale bar represents 10 base substitutions per 100 nt 
positions. The 16S rRNA gene sequences retrieved in this study 
were deposited in Genbank under accession numbers KJ546108- 
KJ546118, KJ588879-KJ588888, KJ588892-KJ588898 (archaeal) 
and KJ588890-KJ588891, KJ588899-KJ588905 (bacterial). 

RESULTS AND DISCUSSION 
FEATURES OF DATASETS 

In this study several metagenomic datasets were examined 
(Table 1). The thalassosaline waters of the salterns in Santa Pola 
and Isla Cristina have a marine origin but they have a different 
source: the Mediterranean Sea and the Atlantic Ocean, respec- 
tively. Besides their high salt concentrations, these environments 
are subject to strong solar irradiation (Rodrlguez-Valera et al.. 



1985; Rodriguez-Valera, 1988). In order to determine how a gra- 
dient of salinity influences microbial communities, we compared 
metagenomes datasets of Santa Pola and Isla Cristina salterns 
and DCM3 and MM5 from marine sites (Figure 1). We expected 
marine derived datasets to form their own branch separate from 
the saltern datasets and that IC21 should be closer to SS19 than 
SS33, because salinity has been identified as the main factor 
determining the distribution of prokaryotic organisms in aquatic 
systems (Lozupone and Knight, 2007; Schapira et al., 2009). 
As expected, the phylogenetic tree showed that marine datasets, 
DCM3 and MM5 clearly differed from the saltern datasets, 
sharing a low number of sequences with them, and the saltern 
datasets were more similar to each other, making it possible to 
recognize the impact of a salinity gradient. Unexpectedly how- 
ever, the community from IC21 was closer in structure to both 
the datasets, SSI 9 and SS33, although the salt concentration in 
IC21 was nearer to the dataset SSI 9, suggesting community com- 
position is affected by local environmental characteristics. Thus, 
we focused this study on the intermediate salinity ponds and 
analyzed the datasets qualitatively and quantitatively to elucidate 
potential causes that might produce the observed differences. 

A well-known adaptive feature for living at high salt concentra- 
tion is the enrichment of acidic amino acids in proteins allowing 
them to properly function at high cytoplasmic salinities (Soppa, 
2006). Therefore, we analyzed the isoelectric point and amino 
acids frequencies of proteins, which would reflect this adaptation. 
Surprisingly, IC21 proteins were more similar to those in SS33 
than to SS19 (Figure 2). The dataset IC21 showed an increase in 
acidic amino acids compared to SSI 9 that therefore indicates a 
greater presence of microorganisms using the "salt-in" strategy 
to osmotically balance their cytoplasm with their environment 
(Oren, 2008, 2013). 

High G-l-C content is often associated with the presence of 
haloarchaea, except for the well-known exception Haloquadratum 
and is a useful predictor of community composition. Analysis 
showed that the G4-C content of IC21 (Figure 3) had a bimodal 
distribution more similar to SS19 than to SS33, and the 
predominant peak in both is at ~65%. This peak is consistent 



— SS37 
-SS33 

— IC21 



-SS19 



-SS13 
-MM5 



-DCM3 



FIGURE 1 I Unrooted neighbor-joining phylogenetic tree based on a 
distance matrix calculated using "bit-scores" between the different 
metagenomes. DCM3: Deep Chlorophyll Maximum, Mediterranean Sea 
(3.8% salinity), MM5: Mar Menor Coastal Lagoon (5% salinity), IC21: Isla 
Cristina saltern (21 % salinity), SS13: Santa Pola saltern (13% salinity), SS19: 
Santa Pola saltern (19% salinity), SS33: Santa Pola saltern (33% salinity), 
and SS37: Santa Pola saltern (37% salinity). 



www.frontiersln.org 



May 2014 | Volume 5 | Article 196 | 3 



Fernandez at al. 



Prokaryotic community of saltern ponds 




SS19 
-IC21 
-SS33 



FIGURE 2 I Isoelectric point (pi) profiles (A) and amino acid frequency (B) of the predicted proteins in the metagenomic datasets SS19, IC21, and 
SS33. The pi was computed for each translated read and is shown as a percentage of the dataset in intervals of bin width 2. 
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FIGURE 3 1 Comparison of G+C% of sequences from metagenomic 
datasets SS19, IC21, and SS33. G+C% was computed for each read and 
the percentage of the dataset in intervals of bin width 10 is shown. 



with the high G+C content associated with most halophilic 
archaea and bacteria described so far (Paul et al., 2008). Ponds 
from Santa Pola saltern (SS19 and SS33) exhibited a low G+C 
peak at 47.9%, which has previously been observed to increase 
at higher salinities (Ghai et al., 2011). This low G+C peak 
comes from Haloquadratum walsbyi (Bolhuis et al., 2006) and 
perhaps from the newly reported nanohaloarchaea (Ghai et al., 
2011). In IC21 the low G+C peak is shifted to a value slightly 
higher than that observed in metagenomic datasets from Santa 
Pola ponds, SS19 and SS33, around 51-52%. This low G+C 
peak might correspond to genera containing representatives of 
halophilic bacterial genera such as Halomonas, Salimicrobium, 
or Salinicoccus. These genera have DNA G+C contents in the 
range of 52.0-74.3, 44.9-51.5, and 46-51.2%, respectively (de la 
Haba et al., 2011). Another possibility is that it corresponds to 
halophilic or halotolerant microorganism not yet described, or 
to AT-rich regions in haloarchaeal genomes (Ram Mohan et al., 
2014) In members of the order Halobacteriales a "minor compo- 
nent" of the DNA (10-30% of the total DNA) with a G+C range 
of 51-59 mol% has been reported (Grant et al., 2001). 

PROKARYOTIC COMMUNITY STRUCTURE 

The taxonomic diversity was analyzed carrying out a search of 
the metagenomic sequences related to the 16S rRNA gene using 



the RDP database and selecting those sequences with a minimum 
length of 100 bp and an identity over 80% for higher taxonomic 
levels (Figure 4) and 95% for the genus level (Table 2). 

In Figure 4 it is observed that the bacterial community 
decreases sharply in IC21 compared to SSI 9, shifting to a largely 
archaeal community. Also identified are members of 14 higher 
taxa in SSI 9 but interestingly we found only six in IC21 and 
five in SS33. The phyla Euryarchaeota and Bacteroidetes and the 
class Gammaproteobacteria are shared by the three datasets, but 
there are numerically more sequences related to the phylum 
Euryarchaeota with increasing salinity and a concurrent decrease 
for the other two. Further, these data show that the phyla biodiver- 
sity in IC21 is more similar to SS33 than to SS19. Ghai et al. (201 1) 
analyzed the changes of the biodiversity along a salinity gradient 
in two metagenomic datasets (with 19 and 37% salts) from Santa 
Pola saltern. The biodiversity detected in the crystallizer pond 
(SS37) is quite similar to that determined for SS33, due to their 
extreme salinities; however in the crystallizer pond only two phyla 
were found, corresponding to Euryarchaeota and Bacteroidetes. 

The simplification of the prokaryotic community at higher 
salinities is also seen at the genus level: the number of sequences 
related to different genera decreases from 69 in SS19 to 28 and 
16 in IC21 and SS33, respectively. The most abundant sequences 
in all three datasets are related to the archaeal genera Halorubrum, 
Haloquadratum, and Natronomonas, recruiting more Halorubrum 
sequences in SS19 and IC21 (12.5 and 65.8%, respectively) and 
Haloquadratum in SS33 (29.5%) (Table 2). The genus Salinibacter 
(belonging to the Bacteroidetes) is the next taxon highly repre- 
sented in Santa Pola datasets with 6.4 and 4.7% of the sequences 
in SS19 and SS33, respectively. However, in IC21 the second pre- 
dominant genus is Psychroflexus, also a member of the phylum 
Bacteroidetes, at 4.6% of the sequences. The six species comprising 
the genus Psychroflexus {P. gondwanense, P. halocasei, P. salinarum, 
P. sediminis, P. torquis, and P. tropicus) have been characterized 
as slightly or moderately halophilic bacteria and most of them 
have been isolated from saline environments (Bowman et al., 
1998; Donachie et al., 2004; Chen et al, 2009; Yoon et al, 
2009; Seller et al, 2012). With respect to genera of the class 
Gammaproteobacteria, in SSI 9 and IC21, Spiribacter is identified 
as one of the most abundant. In spite of the strain "Spiribacter 
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SS19 IC21 SS33 



FIGURE 4 I Representation of high taxonomic levels affiliated to metagenomic rRNA reads. Assigned sequences liave an identity over 80% and a 
minimum lengtli of 100 bp. 



Table 2 | Microbial diversity at genus level in the metagenomic datasets SS19, IC21, and SS33. 



Sample 


SS19 




IC21 




SS33 






Halorubrum 


12.5% 


Halorubrum 


65.8% 


Haloquadratum 


29.5% 




Haloquadratum 


7.9% 


Psychroflexus 


4.6% 


Halorubrum 


23.1% 




Salinibacter 


6.4% 


Natronomonas 


3.2% 


Natronomonas 


5.7% 




Natronomonas 


5.6% 


Haloquadratum 


1.9% 


Salinibacter 


4.7% 




" Spiribacter" 


4.4% 


" Spiribacter" 


1.1% 


Haloplanus 


3.4% 




Haioredivivus* 


3.2% 












Puniceicoccus 


1.9% 












Arhodomonas 


1.8% 












Aquiluna* 


1.5% 












Oceanicola 


1.1% 










Total number of sequences 




1477 




1538 




888 


Total number of genera 




69 




28 




16 



Assigned sequences have a 16S rRNA Identity over 95% and a minimum length of 100 bp. Only those genera with more than 1 % of assigned sequences are shown. 

SS19: Santa Pola saltern (19% salinity), IC2V. Isia Crlstlna saltern (21 % salinity), and SS33: Santa Pola saltern (33% salinity). 

"Candidatus. 



salinus" M19-40 being isolated from Isla Cristina saltern (Leon 
et al., 2013; Lopez-Perez et al., 2013), it was more abundant in 
SS19 compared to IC21 (4.4% vs. 1.1%, respectively). In fact, 
Ghai et al. (201 1) reported a great abundance of sequences belon- 
ging to Gammaproteobacteria related to the genus Alkalilimnicola, 
which later was assigned to "S. salinus" M19-40 (Lopez-Perez 
et al., 2013). Although at the phylum or class level IC21 is more 
similar to SS33 at the genus level the Isla Cristina dataset differs 
in biodiversity and abundance from both Santa Pola datasets, 
demonstrating an overrepresentation of the genus Halorubrum 
and decrease in the sequences (0.4%) related to Salinibacter in 
IC21. Ghai et al. (201 1 ) observed that in Santa Pola datasets (SSI 9 
and SS37) the species Salinibacter ruber appeared as an abundant 
microorganism, but in IC21 it is much lower. 

Traditional culture methods carried out in Santa Pola 
saltern determined that in intermediate salinity ponds there 
were a variety of moderately halophilic microorganisms 
(Rodriguez- Valera et al., 1985). Subsequent studies performed 
by molecular techniques indicated that these prokaryotic 



representatives belong in the groups Gammaproteobacteria, 
Bacteroidetes, and Halobacteriaceae and in a pond with 22% 
total salts 16S rDNA sequences were related to the genera 
Psychroflexus, Halorubrum, and Natronomonas (BenUoch et al., 
2002); similar results were obtained in the dataset IC21. In a 
recent study of a pond from the Exportadora de Sal (ESSA) 
evaporative saltern in Guerrero Negro (Mexico) with 18% total 
salts a Halorubrum-like sequence nearly identical (>99.5% 
similar) to environmental sequences from the Santa Pola saltern 
as well as sequences 97% similar to "S. salinus" were reported 
(Dillon et al, 2013). Recently, Podell et al. (2013) found that the 
microbial composition from Lake Tyrrell (Australia) was corre- 
lated with concentrations of potassium, magnesium, and sulfate, 
but not sodium, chloride, or calcium ions. Sequences related to 
Haloquadratum were positively correlated with potassium, mag- 
nesium, and sulfate ions while sequences related to Halorubrum, 
Haloarcula, Halonotius, Halobaculum, and Salinibacter were 
negatively correlated with them. In addition, H. walsbyi shows a 
higher tolerance to Mg^+ than other halophilic archaea (Bolhuis 
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et al., 2004; Burns et al, 2004). The differences found among Isla 
Cristina and Santa Pola datasets, mainly due to the dominance of 
sequences related to Halorubrum in IC21 and to Haloquadratum 
in SS19 and SS33, might be explained by the difference in the 
ionic composition of both samples as previously reported in 
other hypersaline habitats (Pagaling et al, 2009; Grant et al., 
2011; Boujelben et al, 2012; Podell et al, 2013). 

CONTIGS OF CONCENTRATOR PONDS 

It is expected that contigs assembled from metagenomic reads 
wUl yield genomic fragments derived from the most abundant 
organisms in the sample (Ghai et al., 2011). Assembled contigs 
from the metagenomic datasets SS19, IC21, and SS33 were assem- 
bled and tested against the genome of an abundant organism 
in the datasets, H. walsbyi. In SSI 9, a total of 84 contigs larger 
than 5 kb were assembled, 69 could be assigned to the phylum 
Euryarchaeota and 15 to the class Gammaproteobacteria (Ghai 
et al., 201 1). In IC21 710 contigs with at least 5 kb were obtained; 
618 were assigned to the phylum Euryarchaeota, 85 contigs to 
the class Gammaproteobacteria and 7 contigs to viruses. With 
respect to SS33 a total of 248 contigs were assembled, 247 were 
assigned to Euryarchaeota and 1 to the phylum Bacteroidetes. 
A PCA on the normalized tetranucleotide frequencies of the 
contigs belonging to the most abundant groups was carried out 
(Figures). The contigs from SS19 and IC21 had consistent hits 



to taxa within the Euryarchaeota and Gammaproteobacteria and 
from SS33 the majority were Euryarchaeota. Contigs related to 
the phylum Euryarchaeota are grouped in two different clusters 
with low G+C content, one of them closely related to H. walsbyi, 
comprising 30 and 222 contigs of SS19 and SS33, respectively, and 
the second cluster related to the nanohaloarchaeon "Candidatus 
Haloredivivus" forming 11 contigs of the SS19. Additionally, we 
observed a third cluster including Euryarchaeota of high G+C 
content, comprising contigs from the three datasets that were 
related to genomes of extremely halophilic archaea. On the other 
hand, many contigs related to Euryarchaeota from the three 
datasets did not cluster with contigs related to H. waslbyi nor 
Euryarchaeota of high G+C content. Their G+C content values 
are between the extremely halophilic archaeal reference genomes 
and H. walsbyi; possibly these contigs belong to unknown hyper- 
halophilic archaea. A fourth cluster of Gammaproteobacteria con- 
tigs from SS19 (13 contigs) and IC21 (71 contigs) was closely 
related to the genome of "Spiribacter salinus" M19-40 (Lopez- 
Perez etal, 2013). 

Because of the unique G+C content, the contigs with a 51- 
52% G+C in IC21 were further examined (Figure 3). Contigs 
associated with Gammaproteobacteria and viruses did not have 
this G+C content, ruling them out as contributors. However, 
among the nine contigs related to Euryarchaeota, three could 
clearly be assigned to Halorubrum and also had similar G+C 
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FIGURE 5 I Principal component analysis of tetranucleotide 
frequencies of assembled contigs from SS19, IC21, and SS33 
datasets. Reference genomes are shown as larger circles. The 
following types of contigs are shown: orange: SS19 contigs. Green: 



IC21 contigs. Blue: SS33 contigs, Square: Gammaproteobacterial 
contigs, Triangle: Euryarchaeota contigs. The total number of contigs 
for each cluster, the total length, mean length, and G-l-C% range is 
also indicated. 
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content. These contigs mainly contain genes for ABC trans- 
porters, metallophosphoesterase, multi-sensor signal transduc- 
tion histidine kinase, and hypothetical proteins. 

Further taxonomic analysis was performed on 16S rRNA genes 
from assembled metagenomic reads longer than 500 bp. We found 
eight 16S rRNA sequences from each of the datasets SS33 and 
IC21, and 11 sequences from the dataset SSI 9 that were ana- 
lyzed using BLAST (Figure 6). One assembled sequence from 
each of the datasets grouped within the Haloquadratum cluster 
together with uncultured archaeal sequences from Lake Tyrrell, 
VIC, Australia (with 29% salinity). The 16S rRNA sequences 
from SS19, IC21, and SS33 had high similarity with H. wals- 
byi HBSQOOl of 100, 99.9, and 98.9%, respectively Only one 
sequence from SS33 was found within the Haloplanus cluster 
with a similarity of 98.0% to Haloplanus natans and 97.8% to 
an uncultured archaeon sequence of a saline soil from Jiangsu 
(China). Three sequences from IC21 and one sequence from SS19 
were included into the Halorubrum cluster. One sequence from 
IC21 shared with Hrr. chaoviator a similarity of 99.2%, another 
sequence from IC21 had a similarity of 99.4 and 99.3% with two 
different sequences of uncultured archaea from Aran-Bidgol salt 
lake (salinity over 30%) and other sequence from IC21 shared 
a similarity of 99.5% with a sequence from SS19 and with Hrr. 
orientale. Within the Halohasta cluster one sequence from IC21 
had a similarity of 98.8% to Halohasta litorea. The "Candidatus 
Haloredivivus" cluster was represented by one sequence from 
SSI 9. The rest of the contigs could not be classified into known- 
archaeal genera. The assembled 16S rRNA sequences from the 
three datasets showing the presence of potential microorgan- 
isms related to Haloquadratum (although in IC21 representa- 
tives of this genus do not appear to be very abundant) and 
uncultured archaea in clusters 3 and 5. In IC21 the assembled 
sequences show a high abundance of members of the genus 
Halorubrum. A pattern regarding the community composition 
among the datasets studied appears to be absent, but we do find 
evidence that some sequences from different hypersaline environ- 
ments are related to them. This suggests unknown environmental 
parameters affect community composition, and/or perhaps ran- 
dom dispersal and horizontal gene transfer of key adaptive genes 
play a role. For instance, Parnell et al. (2010) demonstrated that 
genes providing adaptation to their niche, rather than the taxa 
living there, structured halophilic communities in the Great Salt 
Lake. 

Additionally, we found 16S rRNA assembled metagenomic 
reads longer than 500 bp related to the phylum Bacteroidetes: four 
from SS19, two from IC21 and one from SS33 (Figure 7). We used 
BLAST to search in the nr/nt database the sequences showing 
a higher similarity to the 16S rRNA assembled sequences of 
Bacteroidetes. A first cluster related to Psychroflexus was detected 
including one sequence from IC21. This sequence was similar to 
an uncultured bacterium from a Tibetan hypersaline lake (96.4%) 
and Psychroflexus sediminis (95.7%). One sequence from SS19 
was within the Salinibacter cluster, showing a 92.1% similarity to 
Salinibacter ruber. The rest of the sequences from SS19, IC21, and 
SS33 could not be classified into any cluster containing cultured 
microorganisms. The 16S rRNA assembled metagenomic reads of 
Bacteroidetes show that in SSI 9 this taxon is more abundant than 



IC2 1 and SS33, with predominance of sequences related to uncul- 
tured bacteria, while in the dataset IC21 Psychroflexus was the 
genus that recruited more 16S rRNA sequences (Table 2). Neither 
of the 16S rRNA assembled sequences related to Bacteroidetes in 
IC21 are related to any of the assembled sequences of SSI 9 and 
SS33. 

Although at high taxonomic levels IC21 is more similar to 
SS33 than SS19, the fact is that IC21 and SS33 are different 
with respect to the genera and their abundance observed as well 
as to the presence of different uncultured taxa. In saline and 
alkaline lakes in Iran a similar composition of microbial com- 
munities but differing in community structure has been reported 
(Makhdoumi-Kakhki et al, 2012). 

METABOLIC PROFILE 

To analyze the metabolic diversity in the dataset IC21, we 
determined the relative abundance of individual genes involved 
in metabolic pathways of the datasets SS19, SS33, and IC21 
(Supplementary Table 1) by searching predicted metagenomic 
proteins with the UniprotKB database (UniProt Consortium, 
2014) using BLASTX (Camacho et al., 2008). 

Metabolism based on energy from light was queried in our 
data sets. Because our screening protocol removed eukaryo- 
tic microorganisms, only prokaryotic microorganisms were 
considered in this analysis. The photosynthetic reaction centers, 
psbA and pshD genes, were not detected in the three datasets 
analyzed and therefore we can conclude that photosynthesis by 
prokaryotes was absent in these datasets. Instead, genes coding 
different types of rhodopsins were found suggesting that light 
is widely used as an energy source in these conditions, just 
not for carbon fixation. It was observed that bacteriorhodopsins 
and halorhodopsins increased in frequency at higher salini- 
ties. Bacteriorhodopsins associated with Haloquadratum were 
in a greater proportion in Santa Pola saltern datasets while in 
IC21 they were associated with Halobacterium and Halorubrum. 
Halorhodopsins related to Halobacterium were found in high pro- 
portion in all datasets but the number of sequences belonging 
to Haloquadratum is greater in Santa Pola saltern and those 
from Halorubrum are higher in IC21. Sequences related to bac- 
teriorhodopsins and halorhodopsins show clearly the abundance 
of Haloquadratum in Santa Pola SSI 9 and SS37 datasets and 
Halorubrum in Isla Cristina IC2 1 dataset. 

Microorganisms under osmotic stress conditions use different 
survival strategies. Most bacteria maintain cell integrity through 
accumulation of compatible solutes ("salt-out" strategy). 
Sequences related to genes of compatible solutes were mainly 
glutamate synthase, betaine transporters, glycerol kinase, and 
glycerol-3-phosphate dehydrogenase, and at lower frequencies 
were glycerol and glutamate transporters. In IC2 1 lower number 
of sequences related to choline dehydrogenase, glutamate syn- 
thase, and trehalose synthase were observed compared to Santa 
Pola datasets, suggesting a higher synthesis of these compatible 
solutes in Santa Pola saltern. The proportion of sequences 
related to compatible solutes decreased with the salinity except 
for glycerol degradation (glycerol kinase, glycerol-3-phosphate 
dehydrogenase, and dihydroxyacetone kinase), increasing the 
number of sequences in the dataset SS33 compared to SS19 
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Haloquadratum cluster 



.S^Haloquadratum walsbyi HBSC^QO^ [AM1800S8} 
|SS33 uncultured archaeon clone 4099 (KJ588891) 
^IC21 uncultured archaeon clone 1439 (KJ588887) 

Haloquadratum walsbyi C23^ (FR46099) 

- Lake Tyrrell uncultured archaeon clone CBA4087.b1 {JX881 062) 
'-Lake Tyrrell uncultured arctiaeon clone CBA3977.b1 (JX881158) 

-SS19 Ui: 

- Lake Tyrrell uncultured archaeon clone CBA41 23. b1 (JX881 1 79) 
-Death Valley National Park uncultured archaeon clone DV2BarchaeaG01 (AB534147) 

-Dead Sea uncultured archaeon clone DS_ENV_2007_2 {GQ861376) 
-Lake Tyrrell uncultured archaeon clone ABA805-b1 (JX880925) 



jTunis 
-J£l£|Tunis 



n uncultured archaeon clone 186ZA08 {CU467256) 
n uncultured archaeon clone 186ZA05 (CU467259) 



Uncultured haloarchaeon 
cluster 1 



Uncultured haloarchaeon 
cluster 2 



Halogeometricum borinquense PR 3 (AF002984) 
Halosarcina pallida BZ256' (EF055454) 

Guerrero Negro saltern uncultured archaeon clone P11_3_8B (KF234318) 
SS33 uncultured archaeon clone 1515 (KJ588892) 
SS19 uncultured archaeon clone 14449 {KJ546115) 
Australian saltern uncultured archaeon clone Cry7_clone24 (GQ374967) 
Lake Tyrrell uncultured archaeon clone AAA390.b1 (JX880673) 
Xinjiang salt lake uncultured archaeon clone XKL75 (JN14459) 
Cheetham Salt Works uncultured archaeon clone CSW2.24.4 (HQ451075) 
IC21 uncultured archaeon clone 3345 (KJ588882) 
Halobellus davatus TNN1 8' (GQ282620) 

SS33 uncultured archaeon clone 756 (KJ588893) 
Jiangsu saline soil uncultured archaeon clone LA52 {EF533950) 
Hatoplanus natans RE 1 01^ (DQ41 7339) 
Haloplanus aeroganes TBN37^ (GQ282625) 
-Hatoplanus vescus R05-8' (EU931578) 
- Haloferax volcanii HCMB' (AY425724) 




Uncultured haloarchaeon 
cluster 3 



Haloplanus cluster 



- Halopelagius inordinatus R05-2' (EU887284) 

63 Halorubrum chaoviator HALO-G*'^ (AM04B786) 

sj" Halorubrum ezzemoulense 5.1' (DQ1 18426} 

^1' — Halorubnim tebequichense CECT 531 7' (FR870448} 
r-'^^ IC21 uncultured archaeon clone 2064 {KJ588886) 

J — ^ ' Halorubrum distributum 1 m' (D63572) 

' — Halorubrum terrestre 4p' (AB090169) 
'fl ^ Halorubrum arcis AJ201 ' (DQ355793) 

_22 ' Halorubrum lacusprofundi AJCC 49239' (CP001365} 

ea- Aran Bidgo! Salt Lake uncultured archaeon done BC4 (HQ425060) 
J'- Aran Bidgol Salt Lake uncultured archaeon clone BC2 (HQ425055) 
^ IC21 uncultured archaeon clone 7153 {KJ588884) 
89 IC21 uncultured archaeon clone 4188 (KJ588885) 

Halorubrum onentale EJ-52" (AM235786) 

I SS33 uncultured archaeon clone 228 (KJ588894) 

lisian saltern uncultured archaeon clone TSHNAA21 (HQ157630} 
-Aran Bidgol Salt Lake uncultured archaeon clone BA19 (HQ425143) 
I' Tunisian saltern uncultured archaeon clone MHNAA24 (HQ1 57583) 
- Lake Tyrrell uncultured archaeon clone A07_scf00295 (KF673166) 
-Halonolius pleroides 1.15.5 (AY498641) 
- Halohasta litorea R30' (HM1 59607) I 

-IC21 uncultured archaeon clone 1134 (KJ588888) Ha/ohasta cluster 

- Halohasta litchfieldiae tADL' (JF421 969) | 

Halobaculum gomorrense JCM9908' (AB477982) 

Haloarcula marismortui CGMCC1 .1 784' {EF645693) 

—Halosimplex carlsbadense 2-9-1' CAB072815) 

r Guerrero Negro saltern sediment uncultured archaeon clone C10_1G4 {EU570131) 
59- Zhoushan saltern uncultured archaeon clone ZS-35-S2 (KC918829) 




Halorubrum cluster 



Uncultured haloarchaeon 
cluster 4 



JfloT'- Jiani 




Uncultured haloarchaeon 
cluster 5 



Jiangsu saline soil uncultured archaeon clone YA6 {EF533953) 
^ I lOQ — SS33 uncultured archaeon clone 2667 (KJ588895) 

I — Guerrero Negro saltern uncultured archaeon clone P9_2G (KF234314) 

^8 SS19 uncultured archaeon clone 32010 (KJ546116( 

r IC21 uncultured archaeon clone 659 (KJ5B8881) 
™- SS19 uncultured archaeon clone 10790 (KJ546113) 

SS33 uncultured archaeon clone 2473 (KJ588896) 

Tibetan plateau saline lake uncultured archaeon clone CEHLW-A36 {FJ1 55652) 

IC21 uncultured archaeon clone 102 {KJ588879) 
Salton Sea uncultured archaeon clone DSFBPENV12ARC_5B (KC465602) 

!£1 Natronomonas moolapensis CSW8.8.11' [AY498645) 

~ Natronomonas pharaonis DSM 2160" (CR936257) 

Halomarina oriensis KeC-1l' 5AB519798) 
Halalkalicoccus tibetensis DS12 (AF4351 12) 
Natrialba magadii MS3' (X72495) 

(^Guerrero Negro saltern u^^^^^^^^^^ Uncultured haloarchaeon cluster 6 

SS33 uncultured archaeon clone 4270 (KJ588897) ] 

- Halomicrobium mukohataei JCM9738' {EF645690) 
Halorhabdus tjtahensisAK-Z iAF07^8" 



Halobacterium salinanim DSM 3754^(AJ496185) 

ion SS33 uncultured arch 'yv; -c 'K.'S'il^BR 

^ — 



- SSI 9 uncultured arc 



Uncultured haloarchaeon cluster 7 



- Salinibacter ruber DSM 1 3855' (NR_074895) 



Candidatus Haloredivivus sp. G17 (AGNTOOOOOOOO) 
■ Carjdidatus Nanosalina sp J07AB43 (KF673178) 
- Candidakis Nanosalinarum sp. J07AB56 (KF673179) 



"Candidatus Haloredivivus" cluster 



FIGURE 6 I Phylogenetic affiliation of the Euryarchaeaota 168 rRNA 
reads of the datasets SS19, IC21, and SS33. Phylogenetic reconstruction 
was conducted by maximum likeliliood (RAxML) witln near full length (1300 
nt) reference 16S rRNA gene sequences from a manually curated alignment 



and highly variable positions masked. Metagenomic reads were added 
without altering tree topology using maximum parsimony criterion in the ARB 
software package. Bootstrap values are indicated above nodes support 
>50%. The scale bar represents 10 base substitutions per 100 nt positions. 



(number of sequences relative to the total number of metage- 
nomic sequences). This result could indicate a higher presence of 
the algae Dunaliella or a higher primary production of glycerol, 
which is the predominant compatible solute in Dunaliella or 
alternatively, the oxidation of glycerol to dihydroxyacetone, that 
has been studied in the species Salinibacter ruber and is used as 
a growth substrate by H. walsbyi and Haloferax volcanii (Elevi 
Bardavid and Oren, 2008; Ouellette et al., 2013). Therefore, 



glycerol and probably dihydroxyacetone are considered as the 
main carbon source and energy for the heterotrophic community 
in salterns (Borowitzka and Brown, 1974; Borowitzka et al, 1977; 
Ouellette etal, 2013). 

With respect to the nitrogen cycle, it seemed to be simplified 
in the saltern datasets analyzed, with a decrease of the number of 
sequences involved in the reduction of nitrate to nitrite by nitrate 
reductase and nitrite to nitric oxide by nitrite reductase in SS33. 
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Psychroflexus 
cluster 



Tibetan lake unculture bacterium clone SINH1044 (HM127889) 
Tibetan lake unculture bacterium clone SINH432 (HM128002) 

IC21 uncultured bacterium contig clone 1758 (KJSSSSSO) 

Psychraffexus sediminisYlM C238 (EU135715) 
Psychroflexus tropicus LAI (AF513434) 
Psychroflexus saiinarum ISL-14 {EU874390) 
Psychroflexus gondwanensis ACAM 44 (M92278) 
Psychroflexus lorquis ACAM 623' (U85881) 
Psychroflexus halocasel WS 4520 (FR14910) 
Zunongwangia sp. 22II14-10F7 (JQ844757) 
Zunongwangia profunda SM-A87' (DQ855467) 

Lutaonella thermophila CC-MHSW-2' (EU287913) 

Negro saltern uncultured bacterium clone P11_12H_R (KF814617} 
Guerrero Negro saltern uncultured bacterium clone P11_10G_R {KF814610} 

Tunisian saltern uncultured bactenum clone 185ZG06 (CU467419) 
Guerrero Negro saltern uncultured bacterium clone P11_3C (KF234388) 

Israel microbial mat uncultured bacterium clone E4aA01 {DQ103640) 
...(in contig 13242 (KJ588905) 
Gillisia mitskevichiae KCTC 12261 (AY576655) 
Salinimicrobium terrae YIM-C338 (EU 135614) 

Polaribacter dokdonensis DSW-5 (DQ004686) 

Polaribacter franzmannii AlCC 700399 (U14586) 
— Tenacibaculum maritimum R2' (AB078057) 
Tenacibaculum litopenaei B-l (DQ822567) 
Guerrero Negro saltern uncultured bacterium clone P11_7A (KF234385) 
Guerrero Negro saltern uncultured bacterium clone P11_5D (KF814641) 
IC21 uncultured bacterium clone 2812 (KJ588889) 
Israel microbial mat uncultured bacterium E4aF11 (DQ1 03639) 
OwenwBel<sia hongkongensis DSM 17368' {CP003156) 




Uncultured bacteria 
cluster 1 



Uncultured bacteria 
cluster 2 



Alkaliflexus imshenetskiiZ-7Q-\0^ (AJ784993) 

631" Inner Mongolia salt lake uncultured bacterium clone EJ95 (FM210934) 
saL ■ : uncultured bacterium clone 17704 (KJ58890'1) 

^M' -Guerrero Negro saltern bacterium done P12_11E (KF814537) 

d '—Aran Bidgol salt lake uncultured bacterium done OB68 {HQ425210) 
Lj- Guerrero Negro saltern uncultured bacterium clone P11_2C {KF234394) 
''^SS33 uncultured bacterium clone 2054 (KJ588899) 
96| Guerrero Negro saltern uncultured bacterium clone P11_1E (KF814645) 
-SePGuerrero Negro saltern uncultured bacterium clone P11_4E (KF814650) 
' — Aran-Bidgol salt lake uncultured bacterium clone OB95 (HQ425211) 

Cand/datws Amoebophilus asiaticus 6a2 {NR_074574) 

Adhaeribacter aqualicus MBRG1 .5'(AJ626894) 

Gracilimonas tropica CL-CB462' (EF988655) 
Gracilimonas mengyeensisYM J14 [KC561141) 

Balneola vulgaris 1 3 IX/A0 1 /1 64 ' ( AY576749) 

Lake Tyrrel uncultured bacterium done CAB311350411 {JX883442) 
SS19u :]m clone 427 (KJ588901) 



Uncultured bacteria 
cluster 3 



~"Spiribacter salinus" M19- 



Salmibacter ruber DSM 13855' (CP000159) 
Death Valley salt pan uncultured bacterium clone 27-P2 {AB533950) 
11 — Salinibacter iranicus Cb7 (HQ197982) 
'-Salinibacterluteus DGO (HQ197983) 
tO'(CP005963) 



Salinibacter 
cluster 



FIGURE 7 I Phylogenetic affiliation of the Bacteroidetes 16S rRNA reads 
of the datasets SS19, IC21 and SS33. Phylogenetic reconstruction was 
conducted by maximum likelihood (RAxML) with near full length (1300 nt) 
reference 16S rRNA gene sequences from a manually curated alignment and 



highly variable positions masked. Metagenomic reads were added without 
altering tree topology using maximum parsimony criterion in the ARB 
software package. Bootstrap values are indicated above nodes support 
>50%. The scale bar represents 10 base substitutions per 100 nt positions. 



In salterns, sulfate is concentrated along the ponds until its 
saturation and precipitates forming calcium sulfate (gypsum) 
(Landry and Jaccard, 1984). Dissimilatory sulfate reduction has 
been reported until salt concentrations of 24% NaCl (Oren, 
1988). We detected sequences of genes involved in a complete dis- 
similatory sulfate reduction (sulfate adenylyltransferase, adeny- 
lylsulfate kinase, phosphoadenylylsulfate reductase, and sulfite 
reductase) in all datasets studied, except for adenylylsulfate kinase 
in SS33 dataset. The microbial communities of salterns are tra- 
ditionally considered heterotrophic, but chemolithotrophic bac- 
teria can also be abundant and active in extreme conditions. In 
particular, chemolithotrophic sulfur-oxidizing bacteria are able 
to adapt well to hypersaline conditions because the complete 
oxidation of sulfide or thiosulfate to sulfate has high energy 
efficiency (Oren, 1999). In spite of this, we only found some 
sequences of sulfide dehydrogenase that oxidizes sulfide to sulfate 
in SS19. Therefore, an incomplete cycle of sulfate is observed in 
the datasets. In fact, salterns are generally considered eutrophic 
media and so the pressure of selection that favors organisms 
with biosynthetic capabilities with full paths is probably weak 
(Rodriguez- Valera et al., 1981), Additionally, this pathway may be 



carried out by phototrophic sulfur- oxidizing bacteria in anoxic 
sediments, as hydrogen sulfide is considered an important trans- 
porter of electrons between the aerobic and anaerobic habitats 
(Jorgensen, 1982). 

Phosphate regulon (Pho) plays a key role in phosphate home- 
ostasis, products are involved in the transport and use of several 
forms of phosphates (Torriani and Ludtke, 1985; Shinagawa et al., 
1987; Wanner, 1987, 1993). In the datasets studied a scarce num- 
ber of sequences related to genes included in Pho regulon, as 
phoR (environmental phosphate sensor) and phoB (regulon acti- 
vator) were found. However, in the dataset IC2 1 more sequences 
corresponding to the negative regulator protein of Pho regulon, 
PhoU and less for genes involved in the utilization of phospho- 
nate were detected compared to the datasets SS19 and SS33. In 
Santa Pola saltern the total phosphorus concentration increases 
with the salinity (Rodriguez- Valera et al., 1985). However, the 
brines support high Mg^"*" concentrations limiting the availability 
of inorganic phosphate (Bolhuis et al., 2006) and induce the use 
of phosphonate (Fox and Mendz, 2006). A recent study suggests 
the utilization of DNA as a phosphate source (Chimileski et al., 
2014). 
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Overall, the main differences found among Isla Cristina and 
Santa Pola datasets at the metabolic level were a higher number of 
sequences related to genes involved in the synthesis of compatible 
solutes (such as choline dehydrogenase, glutamate synthase, and 
trehalose synthase) and in the utilization of phosphonate in Santa 
Pola datasets with respect to IC2 1 . This is related to the microbial 
strategies of haloadaptation to these extreme environments by the 
different microbial communities present of these habitats as well 
as the ionic composition of the samples. 

CONCLUSIONS 

Santa Pola saltern was built in 1890 over an ancient freshwater 
lake and close to the Mediterranean Sea (Dulau, 1983). By con- 
trast, Isla Cristina saltern was buUt in 1955 over wetlands at the 
marsh of the river Carreras in the village of Isla Cristina, closely 
located to food-processing industries (Moreno et al., 2009). Santa 
Pola saltern is subjected to climatic conditions characterized by 
low annual rainfall and moderate temperatures, with little fluc- 
tuation between summer and winter (Rodriguez-Valera et al., 
1985). However, Isla Cristina saltern is subjected to rainfall sea- 
sons, high solar radiation, and larger temperature fluctuations 
between day and night (Moreno et al, 2009). Although salinity 
has been considered the main factor involved in the structure of 
microbial biodiversity in saline aquatic systems, the data from 
our study strongly suggests that other factors may influence 
the composition of hypersaline aquatic microbial communities 
such as geographic locations (Naor et al., 2012; Zhaxybayeva 
et al., 2013) and environmental characteristics. The environmen- 
tal conditions in Santa Pola saltern are more stable than in Isla 
Cristina saltern, therefore, the prokaryotic community in Santa 
Pola saltern is more stable over time, and in Isla Cristina saltern 
the habitat might be continuously recolonizing. Therefore, the 
differences between these two environments might be due to 
stable or unstable environmental conditions. The IC21 dataset 
looks more similar to SS33 dataset because they are mainly 
composed of representatives of the phylum Euryarchaeota, but 
the reality is that the community structure in IC21 is different 
because the most abundant genus in IC21 is Halorubrum, in con- 
trast to Haloquadratrum, which predominates in SS33 dataset. 
The phylum Bacteroidetes is present in all datasets, but in the 
Santa Pola saltern datasets the most abundant genus of this phy- 
lum is Salinibacter, while in Isla Cristina saltern dataset is the 
genus Psychroflexus. Additionally, in Santa Pola datasets there are 
a higher number of sequences related to genes involved in the 
synthesis of compatible solutes and in the utilization of phos- 
phonate, indicating some differences in the functional activity. 
Despite the results obtained, it is not clear what is causing the 
variations between these salterns; a detailed physico-chemical 
comparative study would be required to elucidate if the microbial 
structure is being influence by abiotic factors or by biogeographic 
situation. 
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